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I. INTRODUCTION 


Electronic Counter Measure (ECM) techniques have gained increasing atten- 
tion in military communications. It is often necessary and important to be able to 
analyze and predict the performance of communication receivers, especially when 
operating under hostile electronic interference. In order to improve the survivabil- 
ity of tactical communication systems operating in Jamming enviornments, spread 
spectrum modulation techniques have been introduced, analyzed and implemented 
in many systems in order to mitigate the effects of intentional electronic inter- 
ference. In many of the analyses carried out to date, the jammer interference 
is typically assumed to be a stationary (or wide sense stationary) random process 
with well-specified statistical characterizations [Ref 1], and the natural interference 
introduced in the communication channel is generally modeled as additive white 
Gaussian noise (AWGN). However, in practice, there are many communication an 
vironments which do not fit the model just described. Consider, for example, a 
communication system which is being jammed by a continuous wave (CW) tone 
near the transmitter operating center frequency, or by a distorted retransmission 
of the transmitter’s own signal. The interference cannot be accurately modeled 
as a stationary random process in either case. Another typical jamming scenario 
involves a jamming signal that may be pulsed between various power levels at a 
rapid rate and is therefore not stationary [Ref 2]. 

The effects of nonstationary interference on the performance of a spread spec- 
trum communication system is not known in general. It is however possible in 
many cases to characterize nonstationary interferers and determine their effect on 


specific communication systems. In fact, if the statistics of the nonstationary 


process — periodic, such a process is said to be cyclostationary (Ref. 3]. The 
statistics and power spectrum of a cyclostationary process are normally computed 
by using phase randomization or time averaging techniques [Ref. 3]. 

In view of the importance of the use of nonconventional forms of generating 
channel interference, the aim of this thesis is to analyze the effects of nonstationary 
narrowband interference on the performance of a spread spectrum communication 
receiver. More specifically, the results of the analysis are applied to the special case 
of a narrowband bi-level pulsed noise jammer, since it is an effective form of ECM 
[Ref 3]. Although certain results concerning the same type of jammer considered 
here can be obtained by other methods [Reference 1 and 2] the results derived are ~ 
based on the application of quasi-stationary assumptions which oftentimes are not 
realistic and eeniiot always be justified in practice. 

In Chapter II of the thesis, methods analogous to those used in the repre- 
sentation of narrowband random processes are used to characterize a stationary 
narrowband nonstationary random process. Such a random process is shown to 
be pacified in terms of its nonstationary quadrature components. The expression 
of the autocorrelation function of the nonstationary narrowband random process 
is then established, and, all the results are used to specify the characterization 
of a bi-level, narrowband pulsed noise jammer which can be obtained by ampli- 
tude modulation of white Gaussian noise which is filtered by a narrowband filter 
centered at some high frequency fp. 

In Chapter III, analysis was carried out in order to characterize the output 
of the spread spectrum receiver consisting of a despreader, a bandpass filter and 
a correlation detector, when a nonstationary narrowband process is applied at the 


input of such a receiver. This characterization is obtained in terms of the mean 


and variance of the Gaussian random variable at the output of the receiver’s deci- 
sion circuitry. 

Chapter IV provides a detailed analysis of the performance of a spread spec- 
trum receiver. The binary phase shift keyed (BPSK) modulation scheme was se- 
lected as the method by which digital information is transmitted and the method 
of spreading the signal bandwidth was chosen to be a direct sequence spread spec- 
trum modulation approach. Noise and jamming were assumed to be sources of 
channel interference, using the models developed in Chapter 3 for their character- 
ization. The receiver performance specification was accomplished by deriving the 
error probability or the bit error rate. The error performance was plotted versus 
various factors that take into account signal power, interference power and band- 
width spreading. The effect of the various parameters on the receiver’s performance 
was analyzed and discussed. 

From the general representation of the narrowband nonstationary random pro- 
cess, other forms of nonstationary interference can be proposed and their effect on 
communication receivers can be derived by similar approaches. Although the de- 
rived performance of the spread spectrm receiver was developed for BPSK signal 
modulation with direct sequence bandwidth spreading, the general results can be 
modified and extended to obtain performance evaluations of receivers in which 


other modulation schemes are employed. 
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Il. NONSTATIONARY RANDOM PROCESSES 


A. NONSTATIONARY NARROWBAND RANDOM PROCESSES 
A random process is said to be narrowband if its power spectral density is zero 
except for a narrow frequency region (W Hz) around a high carrier frequency (f,), 
that is, 
Wie ie ane 21 fous > 1 


If X(t) is a sample function of such a real random process with zero mean, it may 


be expressed as: [Ref. 4] 
X(t) = X, (¢) cos 2a f,t — X, (t) sin 2m ft (2a) 


where the power spectral density of the random processes X,(t) and X,(t) is zero 
except for a narrow frequency region (W Hz) around f = 0. The random processes 
X,(t) and X,(t) are referred to as the quadrature components of X(t). 

Random processes can be generally classified into stationary and nonstation- 
ary processes. If the joint probability density function (p.d.f.) of a process X(t) 
taken at Pa eae times ¢,,t2...¢, is invariant to arbitrary time shifts for any 
integer n, such a stochastic process is said to be stationary [Ref. 5]. A random 
process is nonstationary if it does not have the above properties. A weaker form 
of stationarity is the so-called wide sense stationarity (W.S.S.) involving first and 
second order moments of a random process [Ref. 6]. The representation of Eq. 2.1 
normally implies a W.S.S. random process X(t). 

We propose here to represent a nonstationary narrowband random process in 
a form analogous to that of Eq. 2.1 and to study the statistics of the associated 


quadrature components. This is accomplished by analyzing the output of a 
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narrowband filter when the input is an arbitrary nonstationary process. Before 
attempting this, a convenient method of specifying a narrowband filter is developed 
first. 
1. Lowpass Equivalent of a Narrowband Filter 
Any bandpass filter can be specified in terms of low-pass filter equivalents. 
Let h,(t) be the impulse response of a low-pass filter and H,;(f) be the Fourier 


Transform of h, (¢), that is, 
h,(t) <> Hi (f) = A(fer™ (2.2) 


where H,(f) is not necessarily symmetric about f =0 so that h,(t) is not neces- 
sarily real. Then the Fourier Transform Hz, (f) of a bandpass filter having impulse 


response hg (t) can be specified as follows: [Ref. 7| 


hp (t) <=> He (f) = Ar (f — fo) + He (-(f + fo)) (2.3) 


As shown in Figure 2.1, Hp(f) is realized by summing the frequency up shifted 
(by f,) and the frequency reversed and down shifted (by f,) Fourier Transform 
of the low-pass filter. The negative sign in the term H, (—(f + f.)) is due to the 
frequency reversion of H,(f) described above. 
We can remark that in general h, (t) is a complex quantity. However, since 
the focus is on generating real narrowband random processes, we will require that: 
hg (t) is a real function of t. By inverse Fourier Transformation of Eq. 2.3 and by 


expressing h,(t) in the following complex form, 


hy (t) = hz, (t) + phx: (t) (2.4) 
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H,(f) 





Figure 2.1 Graphic Representation of H,(f) and H,(f) 
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where both h,,(t) and h,,(t¢) are real functions of t, we can obtain the impulse 


response hg (t) of a narrowband filter as [see Appendix A] 
hp (t) = 2Re {hy (t)e*?*4°* } (2.5) 


where Re{(-)} represents the real part of the function in brackets. It is demon- 
strated in Appendix A that h,,(t) must be an even function and h,;(t) must be 
an odd function if h,(t) is to be real. 
2. Output of a Narrowband Filter 

In the case of a stationary random process input, the statistics of the 
output of a narrowband filter is well-known [Ref. 4]. We will investigate the more 
general case where the input random process to a narrowband filter is real but 
nonstationary. Let Y(t) be the output due to a real input signal X(t) applied to 
a narrowband filter having transfer Ametion Hs (f) as described in Section 1 [Fig. 


22\7 Then Y(t) can be expressed as the convolution of X(t) and hg (t) namely 
Y(t) = s ” hp (ta) X(a)de | (2.6) 
It is demonstrated in Appendix B that Y(t) can ¥ a: expressed as 
Y (t) = 2X, (t) cos 2af,t — 2X, (t) sin 27 f,t (2am) 
where 


Xe elie [ [hi, (t — a)cos2a7f,a+h,,;(t —a)sin2rf,a] X(a)da (2.8) 


co 


NG [ [hyi(t — a) cos2rf,a—h,,(t —a)sin2rf,a| X(a)da (2.9) 
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Bandpass Filter 


X(t) Y(t) 





Figure 2.2 Narrowband Filtering of Real Input X(t) 
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There is a certain degree of similarity in the mathematical form of Eq. 2.7 and 
Eq. 2.1. However, the properties of X,(t) and X,(t) given by Eq. 2.8 and Eq. 2.9 
respectively must now be established. 
3. Output Correlation Function for Nonstationary Input 
By definition, the autocorrelation function Ry(t,,t2) of any random pro- 
cess Y(t) is the expected value of the product of Y(t,) and Y (¢,) for arbitrary time 


instants t, and t, {Ref. 5], that is 
Ry (ti, t2) = Ef{y(t:)y(te)} (2.10) 
In order to obtain R,(t,,t2) for the output Y(t) of the filter described in Section 


2, it will be necessary to evaluate the following quantities (since Y(t) is the linear 


combination of X,(t) and X, (t)) 


R, (ti, te) = E{X,(t1)X. (te)} (2.11) 
R,(t,,t2) = E{X, (ti) X, (te)} (2.12) 
R. (ti, te) = E{X, (ti) X, (t2)} (2.13) 
R,.(ti,t2) = Bix. (eee (2.14) 


The four quantities above are complicated functions of t; and tz [See App. Bj. It can 
be observed that in general R,(t,,t2) # R,(ti,t2.) and R,,(t1,t.) 4 —R,.(t1, te); 
unlike the case of stationary narrowband process, both the quadrature compo- 
nents have identical autocorrelation function and the two cross-correlation func- 
tions are null if both X,(t) and X,(t) are real [Ref. 4]. The autocorrelation function 
R, (t1,t2) of the nonstationary narrowband process Y(t) can be formulated as fol- 


lows | See App. B): 


R, (t1,t2) = 4) R(t), t,) cos 2a fot, cos 2aigta tien Cree Sine ee ci ae 
(2.15) 
— Ros (t:,t.) cos2mf,t; sin2myft> —h.. (Cee Sie eos 2 lars 
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It is clear that Ry(t,,¢.) is also a linear combination of R,(t,,t2), R,(ti,t2), Res (ti, t2) 


and ee (t, swe 


B. A MODEL FOR THE NONSTATIONARY INPUT 

The results of the previous sections will be applied now to a specific general 
class of nonstationary random processes which are generated by a system to be 
described now. This is done in order to continue the development for a class of 
nonstationary random processes of interest and to be able to apply the results 
to some practical problems involving jamming of spread spectrum communication 
systems. Let the nonstationary input X(t) described previously be the result of 
mixing a stationary random process W(t) with a deterministic signal q(t) as shown 
in Fig. 2.3. We further assume that the autocorrelation function R,, (t,,t2) of the 


stationary random process W(¢) is 
He inne — E|W (t, )W (t2)] — 5(t. —t,) (2.16) 


That is, W(¢) is assumed to be a white process with unit power spectral density 


(PSD) level. Therefore 
R(t, ,te) = B[X(t1)X(t2)] 
= E (W(t, )a(t.)W (t2)q(t2)] 
= q(t, )a(t2).B [W(t,)W (t2)] 


= q°(t1)5(t2 — t) = 9? (te)6(t2 — t1) 


(2.17) 
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W(t) X(t) 


q(t) 


Figure 2.3 Nonstationary Input X(t) 
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since R,(t,,t2) is zero unless t, = t,. We can now apply the results from Appendix 
B in order to characterize the output of a narrowband filter when its input is a 
nonstationary random process of the type being described in this section. It is 
demonstrated in Appendix C that the autocorrelation terms of Equations 2.11 - 


2.14 simplify significantly when X(t) = q(t)W(t), namely 
R. (t, , £2) = R, (¢, ,¢2) (2.18) 


and 


R.., (t, , t2) = eLearn) — 0 (2.19) 


Finally, the autocorrelation function R, (¢,,t2) for the nonstationary process Y(t), 


namely the output of the narrowband filter, has the following form: [see App. C| 
R, (t,,t2) = 2cos2mf,(t. — t) | hy (t, — a)hy (te — a)q, (a)da (2.20) 
where 


q(t) = 4° (t) (2.21) 


Observe that R,(t,,t.) is not a function of the time difference (t2 — ¢,), so that 


’ Y(t) is (as expected) a nonstationary random process. 


C. A SPECIFIC EXAMPLE 

In the previous sections, we have established a representation in general form 
for a narrowband nonstationary random process generated by applying the mixture 
process q(t)W (¢) at the input of a narrowband filter. As a specific example, it is now 
desired to use these results to model and characterize a pulsed noise jammer that 
ls processed by a narrowband filter. Assuming that we have a bi-level pulsed noise 


jammer, we can model the jamming signal X(t) as the result of mixing white noise 
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W(t) of unit PSD level with a bi-level deterministic periodic signal q(t). Such 
a jamming signal present in the communication channel would cause a typical 
receiver to produce a jamming component at the output of the receiving front end 
narrowband filter. This component is denoted Y(t) as shown in Figure 2.4. 

The autocorrelation function R, (t,,¢2) of the input process to the narrowband 
filter is given by Eq. 2.17 while the corresponding output autocorrelation R, (t, ,t2) 
is given by Eq. 2.20. T, is the period of the periodic signal g(t) which remains at 
the level A for the time interval [o, pT,| and at the low level C for the time interval 
(pT, ,T,],0 < p < 1. [Fig. 2.5] The narrowband filter H, (t) is based on a low-pass 
equivalent which is assumed to have impulse response h, (t) = B- e~®*U(t) where 
U(t) is the unit step function, and B specifies the 3db point of the filter. The 
details of the derivation are presented in Appendix D, where it is demonstrated 


that for the special case under consideration, Eq. 2.20 becomes 





A* B —-B(t,+t2) .2BMin(t,,t3) 
Ry (ts sta) = 2c0s 2 fo(te — t1)) eB tiNe ta) 


B(A? — C?) 
i 
e2BT, . e2BT, [0] — e2BeT, . ——)\ 


eB (tit+ta) panning ats [s| (2223) 
& 


where [z] is the largest integer less than or equal to z 


1 
r= =) + —Min(t,,¢2), (2.23) 
T, 
Ledem 
s=-l1+ 7 Min(t1 , to) (2.24) 


q 


and 


Min(¢,,t2) = minimum of t, and t, 
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Although equation 2.22 is not a simple mathematical form, it is possible to 
derive the average power of the random process Y(t), by time averaging R, (t,t). 


The details are again presented in Appendix D, where it is demonstrated that 
R, (t,t) = B[A?p +c? (1 — p)| (2.25) 


This result can be easily verified from the assumptions made at the beginning of 
this paragraph. Furthermore, it can be seen that for the special case under con- 
sideration, R, (t,t) represents the jammer power at the output of the narrowband 
filter H, (f). 

It has thus been shown that a nonstationary narrowband process can be rep- 
resented by the nonstationary quadrature components in a form analogous to that 
for the representation of wide sense stationary narrowband processes. A general 
expression of the correlation function for a filtered nonstationary input is obtained. 
The results are applied to a special case in order to model a bi-level pulsed noise 
jammer. The effects on a spread spectrum communication receiver due to such a 


pulsed noise jamming signal will be investigated in the next chapter. 
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Ill. EFFECT OF SPREAD SPECTRUM SYSTEM 
ON NONSTATIONARY NOISE 


A. MODEL FOR THE SPREAD SPECTRUM RECEIVER 

In spread spectrum modulation, digital information is transmitted by employ- 
ing a trasmission bandwidth which is much larger than the minimum bandwidth 
required to transmit the digital information via conventional means. In the spread 
spectrin receiver, the demodulation must be accomplished, in part, by correlation. 
of the received signal with a replica of the “spreading” signal used in the transmit- 
ter in order to recover the digital information signal. If the bandwidth is spread by 
direct modulation of a data-modulated carrier with a wideband spreading signal 
of cade the technique is referred to as direct sequence (abbreviated DS) spread 
spectrum modulation. 

The simplest form of DS spread spectrum modulation is DS-BPSK spread 
spectrum transmission where a binary phase shift keyed (BPSK) modulated car- 
rier z spread by a wideband “code signal.” [Ref. 2]. The spreading operation can 
be mathematically represented as a multiplication of the carrier by a function c(t) 
which takes on the values +1 at a digital rate R, periodically. (Fig. 3.1) The trans- 
mitted spread spectrum signal is received together with some type of interference 
and/or Gaussian noise. Demodulation is accomplished in part by remodulating 
with the spreading code c(t). The process is commonly referred to as despreading 
and is a critical function in all spread spectrum systems. 

A simplified model of the DS-BPSK receiver is illustrated in Fig. 3.2. The 


figure shows the received signal r(t) first processed through a front-end bandpass 
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filter centered at the carrier frequency f, and having a bandwidth wide enough 
for the spreading signal (contaminated by thermal noise and perhaps a jammer) 
to be allowed through with negligible distortion. The output of the front-end 
bandpass filter is then mixed with a local synchronized replica of the code spreading 
signal c(t) in order to despread the signal. The despreading operation produces a 
contaminated BPSK signal from which the digital information can be recovered by 


conventional techniques using a filter and a correlator as shown in Fig. 3.2. 


It is well-known that in the absence of any jamming, the perfomance of the 


receiver in Fig. 3.2 in terms of the bit error probability (P,) is given by: [Ref. 2] 
ee 
eS oe 
= 0( 52] (3.1) 


where R, is the bit energy of the received signal, N, is the one-sided power spec- 





tral density level of the AWGN interference and Q(-) is the complementary error 


function defined by 


iE 7 Fat (3.2) 

B. EFFECT OF NONSTATIONARY NOISE ON THE DESPREADER 

AND BIT DETECTOR ON NONSTATIONARY NOISE 

The jamming signal produced by the bi-level pulsed noise jammer may not 
necessarily be narrowband. | However, after the filtering operation performed by 
the receiver’s front-end bandpass filter the jamming signal can be considered to 
be narrowband. We therefore investigate here the effect of the despreader and the 
bit detector system on the input nonstationary narrowband random signal which 
models the jammer. 

The receiver of Figure 3.2 is now assumed to have a narrowband nonstationary 


zero mean input y(t) as illustrated in Figure 3.3. The output of the bandpass filter 
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is denoted by z(t) while the bandpass filter impulse response is denoted h, (t) and 
can be specified in terms of a low-pass filter equivalents denoted h,;(t). The output 
of the receiver due to the input y(t) is denoted Y,. It can be seen that 


Bit) = te ele \ee tert) — [ hy(t —7r)y(r)e(r)dr (3.3) 


oo 


and 


Ts 
Y,; = / z(t) cos 2a fotdt 
0 


which can also be expressed as 


? = | a | / - h;(t —r)y(r)e(r)dr| dt (3.4) 


The jamming of the transmitted signal is always being considered as an action 
of adding a nonstationary Gaussian random process to the signal. As will be seen 
in the analysis to be carried out in Chapter IV, part B, the effect of the jammer on 
the receiver’s performance is completely accounted for by evaluating the variance of 
the output of the bit detector. That is, only the variance of Y,; will be established 
and the other statistical properties of Y; which have no direct influence on the 
systems bit-error probability will be neglected. 

The variance 07 can be shown to be given by an infinite sum of terms, namely 


es 5 ey [Bini (nr, - NPI-APS ~ mR,) 


m= — 0o 


- H, (u)H, (nR, — u) ([S.(f + fo — 4) + S.(f — fo — u)] dudf (3.5) 


where R, is the period of the nonstationary jamming signal and C,, is given by: 


A? p+ c?(1—p) Oo 
C= Ae C2) SIANTP inne 0 
nT p 
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It can be seen from Eq. 3.5 that o% is a function of the duration T, of the 
information bit, the duration T, of the spreading code, chips, and T,, the jammer 
pulsing period. However, practical spread spectrum systems utilize spreading codes 
for which the chip duration is very much smaller than the information data bit 


duration (Ref. 8]. That is 
T, << T, or ieee 


where R, is just T-', the chip rate, and R, is he the bit rate. It is therefore 
possible to define here two classes of jammers, namely fast jammers for which the 
pulse repetition period T, is of the same order of magnitude as T,, so that the 
jammer pulsing rate is much larger than the data rate, and slow jammers for which 
T, is of the same order of magnitude as T,. The effect of the jammer pulsing rate 
on the variance ao; will be examined separately for both cases. 


1. Variance o2 for the Fast Jammer 





Consider first the Fourier Transform P,(f) of the function P,(f) defined 


in Eq. E-2, namely 


P.(f) =| * cos 2m fote-22*!*dt = G(f — fy) + GF +f) (3.6) 


where 


E( f= : Singuehijemian — 


We can see that P,(f) is significant for f = +f> over a restricted frequency band. 
Similarly, P.((f —nR,) is significant for f = +(fp +nR,) over a similarly restricted 
frequency band. The functions P,(f) and P,(f —nR,) are diagramed in Figure 3.4 


from which it is possible to observe that if 
Jo — hy ee Oe 
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or equivalently, if 


2h, << rk, (3.7) 


essentially no spectral overlap of the frequency functions occur. This condition 
except for n = 0 is satisfied for what has been defined as a fast Jammer since 
R, >> R, under fast jamming. 

It is possible to conclude therefore that in the case of a fast pulsed jamming 
environment, the only contribution due to the infinite sum of terms that make up 


the variance o? is the term n = 0 (See Eq. 3.5). Therefore 


o? =C, / ” Hy (f) Hs (-1)P.(-S) 
ae (3.8) 
/ Hy (u) Hz (—u) [Sef + fo — u) + S.(f — fo — u)] dudf 


and substituting in Eq. 3.8 the corresponding mathematical description of the 


functions H;(-),H,(-),P.(-) and S,(-) as given in Eqs. E.3, G.6, G.8 and G.17, 








results in 
Cot 2 aay 1+e Te 
cee OO _ _ p-BT.\ _ 
v7 = 3 {1 i: ( RF az) (1 — 7? **) Bae 
B? B(B? + 8?) 
-bT, 
+ (2e°*7* — 1) Ry mre eypt 9) 


2T,T.(B +b)? | 2bT,T.(B —b) 


Be-2T< e- *T- Ben BT eeTe | 
The details of the derivation of this result are presented in Appendix G. 
2. Variance o? for the Slow Jammer 
In the case of a slow pulsed jamming environment, the condition R, >> 


R, is not satisfied, therefore a certain degree of overlapping between the spectral 


factors that make up the expression for 0} will occur. While the infinite sum of 
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Eq. 3.5 can be approximated accurately by a finite sum, the number nj, of terms 


needed in the approximation will be such that 
Mo ses >: 2h, 


or equivalently 


2R, 
>> — 3.10 
ho R, ( ) 


In this case, the result for 0} becomes 


a= > Of” HSH (ER, —NP.(-DP.US - FR) 


k=— no 
[Ba (w)He (ER, — w) [S(f + fo — uw) + S.(F - fo — u)] dud 
_ (3.11) 
It is also clear that as k increases, |C,,| decreases so that the sum involved 
can merhaps be further truncated. The computation of 07 in this case becomes 
quite complicated, however using numerical integration, 07 can be obtained quite 


accurately. 
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IV. PERFORMANCE OF THE SPREAD SPECTRUM 
COMMUNICATION SYSTEM 


A. GENERAL OUTLINE 

In Chapter III, a complete derivation for 04, the variance of the Gaussian 
random variable that quantifies the power contribution of a bi-level, nonstationary 
narrowband jamming signal at the output of a DS-BPSK Spread Spectrum Mod- 
ulation receiver, has been undertaken. A diagram of the spread spectrum receiver 
is shown in Fig. 3.1. The receiver structure is optimized based on principles of 
Maximum A Posteriori (MAP) detection for binary modulation schemes. The re- 
ceived signal and noise is first processed by the front-end bandpass filter which has 
a bandwidth just wide enough to pass the spread signal and the in-band noise. The 
narrowband signal is then processed by the rere which consists of a multi- 
plier mixing the received signal with a synchronozed local replica of the spreading 
code signal. A second bandpass filter having a bandwidth proportional to the in- 
formation bandwidth follows the despreader, so as to further eliminate the noise 
power being presented to the remainder of the receiver. The data signal is then 
demodulated by a correlator and finally, a decision is produced every T, seconds 
which is translated into O or 1 data bits. This generates a MAP estimate of the 
digital data at the output of the receiver. 

The performance of such a receiver is usually characterized by the probability 
of error, which is, in essence, the probability of an incorrect decision made at the 
receiver. The probability of error can be further categorized into symbol error and 


bit error probability. However, for the case of Binary Phase Shift Keyed (BPSK) 
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modulation, the two error probabilities are equivalent and provide the basis for 
specifying the receiver’s performance. 

In the next section, we investigate the performance for the case in which a 
BPSK modulated signal is spread by direct sequence methods, and the signal is 
received in AWGN interference as well as the bi-level pulsed noise jammer described 


and analyzed in previous chapters. 


B. RECEIVER PERFORMANCE FOR BPSK MODULATION 
The structure of a spread spectrum correlation reciever is shown once again 


in Fig. 4.1. In BPSK modulation, the two transmitted signals are [Ref. 4] 


si (t) = \/ Fe" cos 2a ft O<t<T, and 
b 
2E 
S(t) = —4/ =~ cos 2m fat Oat < T, 
b a 


where £, 1s the energy of each signal and J, is the duration of each signal. Since 


(4.1) 


each signal represents a single information bit, E, is the energy per bit and T, lig 
the bit rate. 

At the output of the second bandpass filter, the signal y(t) can be regarded as 
the summation of the transmitted signal s;(t), narrowband noise N;,(t), and the 
jamming signal y,;(t) convolved with the impulse response h;(t) of the bandpass 
filter. That is, 


y(t) = s (t) + Men (t) + yy We(| Aw (4.2) 


4 = 1] or 2 


Hence, 


Ty 
ve} y(t) cos 2m f,tdt 
0 


_ | . S;(t) *« h(t) cos 2a f,tdt + T, Nin (t)C(t) *h,(t)cos2af,tdt (4.3) 


0 


+ 


/ yz (t)c(t) « h,(t) cos 2 f,tdt 
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or simply, 


ia = S; = Ne + Y; (4.4) 


where S;, N., and Y,; represents the three integrals in Equation 4.3 respectively. 
It is clear that Y is in fact a random variable whose statistics can be evaluated 
once certain properties of S;,N;, and Yj, are established. 

We first assume that the noise contribution at the front end of the receiver is 
zero mean, white Gaussian noise with vo as its two sided power spectral density 


level [Fig. 4.2]. The autocorrelation function for the noise, ednoted Ry (t,r) is 
N 
lt) amok! — 7) (4.5) 


where 6(t) is the Dirac Delta function. Based on these assumptions it is shown in 
Appendix H that N,, is also a zero mean, Gaussian random variable with variance 
o?, specified in the sequel. Due to the fact that the receiver is linear, the noise 
contribution at the output retains its Gaussian statistics. The expression for o7, 


is from App. H given by (Eq. H.18). 


N,T,k 
—— = (4.6) 
where 
k = 0.903 
and 


c A ‘ : 
= processing gain. 


a 
Ry 


It is shown that S, and S, are given by [Ref. 2] 
s, = aa $= (4.7) 
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Figure 4.2 Power Spectral Density of the Thermal Noise 
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Furthermore, Y, is the quantity related to the jammer and its variance 0% is given 
by (see Eq. H.4) 
Cri, 


27 = —J(G,. 
Oy; 8 ( ) 


where J(G,) is defined by Eq. 4.6. To further simplify the notation, let X = 
Nin + Y;7, so that 


e— oF, + C; (4.7) 


since the noise and jamming are zero mean uncorrelated processes, and 


then X is zero mean Gaussian random variable with pdf given by [Ref. 2] 





f.(X) = eae (4.9); 


For the BPSK modulation scheme that is being considered here, the decision 
regions of the signal space diagram can be illustrated as shown in Fig. 4.3 [Ref 4]: 
From equations 4.8 and 4.9, it can be seen that the probability of receiving the 


signal Y given that a signal S; was actually transmitted is: 
Pr{r = Y/S;(t) was transmitted} = f,,,,(Y/S:;) = f2(Y — S;) (4.10) 


The probability of bit error rate is then given by: [Ref. 2] 


ET, 
n= 9( ae 





or 


P, = o( sare) (4.11) 
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Figure 4.3 Signal Space Diagram for BKSP Modulation 
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where Q(-) is the complementary error function defined earlier in Chapter III. 


Substituting 07, and of} from Equations H.4 and H.18 results in 


EAT, 
oe) 
N f k ot CaF J(G.) 
or 
Ey 
P, — o( ae (Cad(Gey CoJ(Ge) (4.12) 
wT WF 4 


C. ANALYSIS OF ERROR PERFORMANCE 

As derived in the previous section, Eq. 4.16 yields the probability of bit error 
P, for a DS-BPSK spread spectrum receiver under AWGN noise interference and 
narrowband bi-level fast pulsed noise jamming. Clearly P, is a function of the 
power spectral density level N, of the AWGN, the average power of the jammer 
given by BC, and also, the processing gain G,. It is clear from Eq. 4.11 that 
the receiver bit error probability will increase if the noise and/or jammer power 
increases as the Q(-) function is monotonically decreasing. 

In order to analyze the quantity P, in more detail, we define the signal to noise 
ratio (SNR) to be the ratio of E, to N, and the jammer power to signal power ratio 


(JSR) to be the ratio of average jammer power P, to the average signal power P 


[Ref. 4]. Hence, 


EB, P,; 
SNR = — wR Alki 
N. Pp (4.17) 
where 
i Cre 2H. (4.18) 
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so that Eq. 4.16 can be reexpressed as 


ous ; (4.20) 


P= ( o903 4 (JSRIISNR) 
* 8nG. 

Based on Eq. 4.20, families of curves of P, versus SNR for fixed JSR values can be 
plotted and analyzed. However, first plots of P, versus JSR for different values of 
G, under various SNR conditions are shown. Logarithm scales are used for P, due 
to the large variation in this parameter. From Figure 4.4a to d, one can observe 
that in general, under intense jamming (i.e., JSR large), the receiver probability of 
error is also large, thus implying that a significant number of delivered data bits 
are in error. Under a fixed JSR and fixed G, values, it can be seen that if the SNR 


is large, P, can be made significantly small. For example in Fig. 4.4b, it is shown 


that when JSR = 15dB and G, =100; 
Paat SNE — 4) = 10r.e 


and 


Plat sik — 20) 10a. 


However, when the jamming is very strong (for JSR > 30 dB) P, tends to a limit of 
approximately 10~* regardless of the value of all other factors. It can be observed 
from Eq. 4.16 that if JSR increases indefinitely, the limit of P, equal to 0.5 is 
reached. 

From Fig. 4.4a to d, we can also conclude that for any fixed values of JSR 
and SNR, an increasing processing gain G, implies a decreasing value of P,. The 
phenomenon is more prominent in a high SNR environment than in a low SNR 
environment. There is no significant improvement in P, for SNR in the range from 


0 dB to 8 dB however a large improvement in P, is obtained when SNR exceeds 
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Figure 4.4 Pb versus JSR for (a) Ge 
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12 dB. The immediate conclusion that we can draw from this is that by increasing 
the chip rate (R,) of the spreading code, resulting in a higher processing gain, 
the perfromance of the receiver can be improved under the jamming environment 
analyzed, however this counter measure loses effectiveness when powerful jamming 
is present. 

To better illustrate the effect of G. on P,, another family of curves is plotted. 
In Fig. 4.5a to d the variation of P, as a function of JSR is shown for different 
values of G, and SNR. It is clearly shown that all the curves converge toward P, 
equal to 0.5 as JSR increases. From Fig. 4.5a we can see that P, (at JSR = 20 


1.5 


dB) reduces from approximately 10-°-’ to 107 *-° when G, increases from 2 to 400 
under weak signal (SNR = 4 dB) conditions, while P, reduces from approximately 
10-°-” to 107° under the same conditions when SNR = 20 dB (see Fig. 4.5d). 
At JSR > 35 dB, no significant gain can be obtained regardless of signal strength. 
As a result of this, we can conclude that by increasing the chip rate to counter the 
narrowband bi-level fast pulsed noise jammer is not effective under strong jamming 
conditions, however, a high chip rate can reduce the probability of bit error of the 
jamming is not too severe. 

It is also possible to represent the error performance curves by plotting log,, P, 
versus SJR, the signal power to jammer power ratio under various SNR conditions. 
Fig. 4.6 a to d show how P, varies with SJR. The curves can be interpreted 
in essentially the same way as the previously mentioned graphs. By comparing 
Fig. 4.6 a to d with the results obtained in Reference 2, we observe that the 
performance of the DS-BPSK spreak spectrum peciever degrades under a jamming 
environment in a manner similar to the case of barrage noise jamming or partial 
band noise jamming environment presented in the reference. Furthermore, the 


receiver performance is well below the boundary of the performance for the 
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(b) Ge 
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worst case pulsed noise jamming derived in Reference 2 based on a quasi-stationary 
analysis method. Except for the weak signal condition case, the receiver perfor- 
mance approaches the above mentioned boundary. In general, by increasing the 


processing gain, the performance of the receiver improves significantly. 
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V. CONCLUSIONS 


In this thesis it has been demonstrated that any nonstationary random process 
that is bandpass filtered can be represented in a mathematical form analogous to 
that used to represent narrowband wide sense stationary random processes. The 
quadrature components of the narrowband nonstationary process are themselves 
nonstationary low pass random processes. In general, the autocorrelation function 
of the two quadrature components are not identical and the cross-correlation func- 
tion of the quadrature components is not necessarily an odd function. However, 
when the input nonstationary random process is real, the correlation functions of 
the quadrature components of the narrowband output process have properties sim- 
ilar to those associated with narrowband wide sense stationary random processes 

The general representation of the nonstationary narrowband random process 
has been applied as a model for a specific class of electronic einitit measure 
signals. Specifically a bi-level pulsed noise jammer, in which the power of the 
jammer periodically pulses between two levels is modeled using the previously 
mentioned nonstationary narrowband random process representation and its effect 
on a DS-BPSK receiver analyzed. This type of jammer was analyzed because it is 
quite effective in the sense that it produces significant performance degradations 
on the receiver. In practical cases, the jammer pulses between two levels at a rate 
faster than that of the data bit rate. The quasi-stationary analysis method for 
recelver performance is not valid. Therefore, a methodology has been developed 
that allows rather precise evaluation of receiver performance in terms of signal, 


noise, and jamming powers as well as spread spectrum processing gain. 


ol 


The effect of the pulsed noise jammer on the spread spectrum binary phase 
shift keyed receiver was examined in detail and it is shown that a strong transmitted 
signal is more effective at mitigating the effect of the jammer than increasing the 
chip rate of the spreading code. The pulsed noise jamming can be shown to be 
quite effective under various signal power and processing gain values, as rather 
large bit error rates result when powerful jamming is present. 

The thesis develops a methodology that can be used to evaluate the effect of 
a nonstationary narrowband random process interference on communication chan- 
nels. There are two benefits from this study. First, in terms of ECM performance 
it makes it possible to evaluate the effectiveness of such jammers, and second, in 
terms of electronic counter counter measures (ECCM), the performance of jam re- 
sistant receivers can be evaluated more accurately. Although the example chosen 
in the study involves BPSK modulation spread by direct sequence methods, results 
for other forms of nonstationary jammings on different types of spread spectrum 


receivers can also be obtained from the above results with appropriate modification. 
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APPENDIX A 
REPRESENTATION OF NARROWBAND FILTERS 


Let 
h(t) 2+ Hi (f) = Ae (4.1) 
be the specification of a low-pass filter (LPF), and let 


hz (t) ——> Hy (f) = Hi (f — fo) + Hi (-(f + fo) (A.2) 


where f, is assumed to be much higher than the cutoff frequency of the LPF. 
Thus, Hz, (f) is the transfer function of a narrowband filter centered at f = +f,. 


By inverse Fourier Transformation, we have: 


hp (é) =| - Hy (f)e??""* df 


=f mu-pemra+ | m(-t-pyerrag 49) 
= hr (a + hr (—t)e~??* fot 
If we express | 
h, (t) = hz. (t) + phic (t) 
since h, (t) may be complex, then 
hz (t) = lor (t) + hr, (—t)] COS Ziehet aaa or (—t) = hr; (t)| sin 2a fot+ 
y {hz ;(t) + hy; (—t)] COs Zon fee + oon (—t) = H,- (—t)] sin Diepat (A.4) 


Since hg (t) will be assumed to be real, the imaginary part in Eq. A.4 must vanish, 


so that 
hy: (t) = —hzi(-—t) =>  h,; is an odd function. (A.5) 
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and 


hy-(t) =hr-(-t) => At, is an even function. 


Observe also that 
h,(—t) = Az, (—t) + hii (—t) = hi, (t) — shi: (t) = h(t) 


Therefore 
hg (t) = 2h,, (t) cos 2af,t — 2h; (t) sin2a fot 


=P (t)e?2" fot i [Ay (t)e2?™ fot) 
and finally 
hg (t) = 2Re {hr Qerwe 


o4 


(A.6) 


(A.7) 


(A.8) 


(A.9) 


APPENDIX B 
OUTPUT OF A NARROWBAND FILTER 


The output y(t) of a narrowband filter due to an input z(t) is 
a) = ™ hp (t — a)z(a)da 
a - [2h,,(t — a) cos2af,(t — a) — 2h,;(t — a) sin 2rf, (t — a)] r(a)da 
~_ | (B.1) 
where the representation derived in App. A has been used to obtain Eq. B.1. By 
expanding cos 27f,(t — a) and sin27f,(t — a) it is possible to obtain 
y(t) = 2X,(t) cos 2af,t — 2X, (t) sin 27 ft Dy eZ) 
where 
XC) a [h_-(t — a) cos2a7f,a+h_,;(t — a) sin2rf,a] X(a)da (373) 
and 


rane) — ae [rz (t — a) cos2mf,a—h,,(t — a)sin2af,a] X(a)da (pae4) 


co 


so that a standard narrowband-type representation is possible for the output y(t). 
If X(t) is a random process that is not necessarily stationary, X,(t) and X,(t) are 
themselves random processes with autocorrelation functions 
Ratz) = &{X.(t:)X.(t2)} 
~ i 1 ; * —j2rfoa 
al [ 5 [Pes (th — ax)e?™22* + hy (tr — a)e7?7/°*| 
1 
2 
mie. [> 
= al / Re {hy (t; — a)h, (te — peitrtola+8) | R, (a, B)dadB+ 


[hz (t2 — B)e??"/°? + ht (t2 — B)e~??*/°?| Rz(a, B)dadB 
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sf [ Re {hi (t, a) hz (to — pei?) | R, (a, B)dadp (B.5) 


and 
R, (t, ito) = E {X, (t:) X, (t2)} 
= “ats I. Re {hs (tr — aha (ta — Be s?r/("*") } Ry (a, B)dadp+ 
of [Ref Re{ hy (tr — a)hy (te — pye**!(2-#) | R, (a, 8) dadp 


(B.6) 


The cross correlation functions are given by 


[> 


R,, = E{X, (t:)X, (tz)} 


sf i Im { fhe (t — a)hy (tz — B)e??*Fe'e- 9) \ R, (a, B)dadB+ 


off. Im | {hi (t, — a)hyz (te — pe #2049) | R, (a, B)dadg (B.7) 
and 
Rett) =e eaten ei) 


af [- Im {hi (ti — ahs (te — B)eW274(*-9) } R, (a, B)dad p+ 


—5 ff tm fii te — adhe (te — Bern”) } R, (a, 6) dade 
(B.8) 
Observe that in general X,(t) and X,(t) normally labeled the quadrature compo- 
nents of the narrowband process y(t) are not stationary in general. Nevertheless, 
the derived autocorrelations and cross correlations can be used to obtain the au- 


tocorrelation of the random process of y(t), namely 
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Ry (t, ste) = E {y (ts )y(te)} 


= Bax, (t,) cos 2af,t, — X,(t,) sin27f,¢, | 
-2|X.(t2) cos 2a f.te — X,(t2)sin 2n fatal} 


=A R. (t, ,t2) cos2mf,t, cos 27f,t. + R,(t,,t2.) sin2arf,t, sin2af,tet+ 
— R,,(t1,t2) cos 27f,t, sin 2x f, te 


— R,,(t,,¢.) sin 27f,t, cos an fats (B.9) 
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APPENDIX C 
NARROWBAND PROCESS AUTOCORRELATION FUNCTION EXAMPLE 


The autocorrelation functions determined in Appendix B are now used to 
treat a specific case, namely when the input process X(t) to the narrowband filter 
is generated by the system of Fig. 2.3. The resulting output y(t) is a narrowband 
process bearing quadrature components whose autocorrelation function is given by 


(see Eq. 2.17) 


R,(ti,t2) = s/ Re {hy (t; — a)hyz (t2 — Ge) ei alae nas (a)da+ 
a (C.1) 
o/. Re {hi (t; — a)h, (te — a)}¢, (a)da 


Given that h, (-) represents a low-pass system and q, (a) would normally vary much 
more slowly than the frequency 2/,, it can be seen that integrals involving a 2f, 


frequency are negligible resulting in further simplifications. Thus 


Re (ty ste) = 3 / ~ Re tht (a) eat; See ae (C.2) 


oo 


and similarly 


atec i sf. Re {hi (¢, — a)hy, (t. — a)}9¢,(a)da = R, (t, , te) (c.3) 


oO 


The crosscorrelations are given by 


oe tl sf Im {hy (t; — a)hy (te — a)}q,(a)da = R,.(t,,t2). (C.4) 


oo 
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For the special case where h, (t) is real, R,,(t,,¢.) and R,,(t,,t.) are both equal 


to zero, so that from Eq. B.9 and Eqs. C.2 and C.3, R,(t;,t2) becomes 


R, (t1,t2) = 4B, (t,,t2) [cos 2af,t, cos2af,tz +sin2af,t, sin27f,te | 
(C.5) 
= AR. (t, ; te) cos 27 fo (t2 = t,) 


and finally 


R, (t, ,t2) = 2cos2rf,(t, — ts) | h,(t, — a)h, (t. — «)q, (a)da (C.6) 


— co 
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APPENDIX D 


AUTOCORRELATION FUNCTION AND AVERAGE POWER AT THE OUT- 
PUT OF A NARROWBAND FILTER DUE TO A BI-LEVEL PULSED JAMMER 


When q(t) is a bi-level signal as shown in Fig. 2.5, it is clear that 


q, (t) = A* — {a ou) s [U(é - kT, = pe) =z Ut - kT, 1} (D.1) 


k==- 00 


where U(t) is the unit step function. Let 
I(t) = | h,(t, — a)h, (t2 — a)q,(a)da 


This integral must be determined in order to evaluate R,(t,,t2) as given by Eq. 
C.6. Substituting g,(t) into I(t) and utilizing the one pole filter specification for 


h,(t) having 3 dB cutoff B rps, we get 


A’B , 
I(t) = 5 e7 Bltitts) 2B Min(ti,t2) — B?(A? _ Cer 2 a) 
 f e?F2U(t, — a)U(t. — a) Se U(a—kT, — pT,)da (D.2) 
ao k=- 00 
-| e°F*U(t, — a)U(t. — a) yy U(a—kT, —T,)da] 
a k=-— oo 


The remaining integrals in Eq. D.2 are: 


oo Min(t,,t3) 1 [9] 
2Ba — 2BMin(t,,t3) 2B(k+1)T 
eo" “da = — le —e 
Z i ee 


k=- oo k—p)T, k=- oo 
D3 
if (k+1)T, < Min(é,,t,) ‘alt 


= 0 otherwise 
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and 


co Min(t,,t3) 1 [3] 
| 2 se anette a Pa(eer Te | 


Mein VIR A 1): 2B es (D.4) 
if (k+1)T, < Min(t,,t2) 
= 0 otherwise 
where [z] = largest integer less than or equal to z and 
_ Le 
r=—e+—Min(t,,t,) s=1+—Min(t,,t.) (D.5) 
le I, 
The difference of the two integrals given by Eq. D.3 and D.4 produces 
1 [r] [>| 
= » Bye teste) a eaterelts| = S- eer — e2B(kt sais 
k=-— © k=>-©0 
(D.6) 
since 0 < p <1, we have rr > s, so that 
[r] [>| 
3 e2BMin(ti,t2) _ ys e2B Min(ti,ts) — e2BMin(ti,t2) ([r] = (s]] (D.7) 
k=—-— 00 k=>- 00 
Furthermore 
[s] [r] 52 BT 2BT,(s| 2B pT 2BT,(r| 
2B(k+1)T, _ AB (RY pili. Cale eae ee © eee 
ye c De - ; a 1 — e-28T, (D.8) 
=- 0 k= — 00 
Combining the above results yields (see Eq. C.6) 
R, (t1,t2) = 2cos 2mf,t - I(t) 
2 
= 2cos 27 f 1 <* Fe- Blta tts) g2BMin(ta.t2) _ 
ane a2 
B(A? — c?) (D.9) 


5 e7 Bltitts) e2erententnd (o] —(] 


Sy ELS (Cd oe ne 


1 — e-253T, 


+ 
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Since X(t) represents a bi-level pulsed jammer, the power due to the jammer 
at the output of the narrowband filter denoted by P, can be obtained by letting 


t, = t, =t in the expression of R,(t,,t2) and time averaging the result. Thus 


2 v2 
R, (t,t) = A? B — B(A? — C?)({r] — [s]) — oa e2BTs 9-28 (t-Tel4]) 


— e-2BT, 
— e2B pTe p— 2B (t-Te[r]) 

(D.10) 
which must now be time averaged in order to obtain P,, the average jammer power. 
With the aid of Diagram D.1, D.2, and D.3 we note that all the three functions 
{(r] — [s]},{t —T,|s]} and {t — T,|r]} are functions of t, and. their time average 


values can be calculated as follows: 





(I= BSD = tim, (elt +(e 1901 - 0) )) 








: 1 
= lim mal +(1-—p)n]=1-p (Dae) 
e~ 2B(t-T,[S]) lim : = »~ 2B (t+t,) aa (1 -2BT (D 12) 
ir ne JU - ‘ 
eee 2BT, ae 


SS ; il eT, T 9 
e 2B -T ID) = lim marl | cain | creer dy 
ae OO (n + p)T, 0 0 


pn DPD (eee 2 a 
= EE (a 
: q 


Therefore, the final expression for P; is 


2 772 _ 2-2 asd 
Py = APB B(A? ~ 07)(1 — p) - FASE) [Pe _ te 


= B[A’?p+C7(1—- > )| 
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(1+0Tq) (2479) 


(sl=[=—1] 





7q Gere (2+prq) t 


Figure D.1 Function [r]—[s] 


63 


z=t 
z=Tq[r] 
37 
2Tq 
Tq 
(14°79) (2+PTq) i 


(2+97q) 





Figure D.2 Function (t-—Tg[r]) 
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Figure D.3 Function (t-—Ta[s]) 
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APPENDIX E 
VARIANCE OF Y; 


Since y(t) is Gaussian and zero mean random process, Y; is a Gaussian, zero 
mean random variable due to the fact that linear operations only are applied to 


y(t). [Ref 2] The variance of Y,, denoted by o4 can therefore be expressed as 


follows: ee : 
= By [ / 2z(t)z(v) cos 27 fot cos 2n fovdtav 
" 0 0 
ter: (E.1) 
= | | R,(t,v) cos 2m fot cos 2m foudtdu 
0 Jo 
where 
Ea | / hz (t= 2)fer (oR foo eles (E.2) 


and R,(r —-y) is defined as the autocorrelation function of the spreading code. For 
DS-BPSK spreading modulation, R,(6) is well known, namely 


6 * 
R(8)={1—% eee 
0 "| Olezaae 


where JT‘ is the digital rate of the spreading code and is called the chip rate. 
The power spectrum of the spreading code S,(f) is the Fourier transform of R, (6), 


namely [Ref. 2]. 





Sah i este a (ee (£3) 
where 
; sin TZ 
sinc) — = 
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— 1, Zo & 


| 5 (f)=sine “(fT ) 





Figure E.1 Autocorrelation Function and Power Spectral 
Density of DS-—BPSK Spreading Code 
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A diagram of R,(6) and S,(f) is shown in Fig. E.1. 


A new signal P,(t) is defined as follows: 


P,(t) = ao “Oa 7, 
7 - otherwise 


So that equation E.1 can be rewritten as: 


_ / - [ i R, (t, v) P, (t)P, (v)dtd 


(E.4) 


It is possible to express R, (t, v) as the double inverse Fourier transform of Sz (f,v), 


that is 


a) — | / Ss? (fi vye 27 S94) day 


~= p [ R, (t, v)e7 2278") dedy 
Let P.(f) be the Fourier Transform of P,(t), so that Eq. E.2 becomes 
a= [- [- Sz (f,v) )f- frm (te 2"! Pp lemsetidtdudidy 
=| f se(tvyPe(—Dpe(v)afav 


Using equation E.2 to evaluate equation E.6 results in 


y= f- [2 (7,9) R.(7 -¥ )f- fs Kleen 


h,(v — -y)e?*""”" dtdudrd-- 


= H,(f)H1(-v) [ [ R,(7,7)R.(r — y)drdy 


/ hy (t— re"?! dt = Hy (f)e-227!7 


where 


since 


68 


(E.5) 


(E.6) 


(E.7) 


(E.8) 


(E.9) 


and H,(fo) is the Fourier Transform of h;(t). Therefore 


=[ [mine (rn. 


(E.10) 
[- [ (7 — ye 2774" 63?" drdydf dv 


The double integration involving the variables 7 and 4 in equation E.10 can be eval- 
uated for the special case where y(t) is the nonstationary random process described 


previously. 


In order to evaluate o*, the double integral D(r,-y) defined as must first be 


evaluated, where 


D(r,1) =f [® (r,q)enF?"t ef?" Mdrdy = (E11) 


Recall from equation 2.20 that 


R, (7,7) = 2.cos 2a fy (r — y)n (7,7) (#.12) 


where 


Fa t= / (ee =a) a: a). (E.13) 


oo 


The double Fourier transform pair n(7,-y) and H(u, w) can hence be expressed as: 


(u, w) =f [x (Geyleme tee Cardy (E.14) 


so that 


n(t,) = / | icreem te © ere! dudw (E15) 
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Therefore, 


vena) = [7 [af [7 [7 mtanjerre ser aude] 


R,(r — y) cos 2m fo(r — ye" 77777 €??*" drdy 


=f- o H(u,w) f [- ila eet (ee 


e~F2e(f—u)t e-52"(w—¥)9 de dadudw 


=f ns H (uw) f IS.(f —ut+fo)+S.(f —u—fo)] 


e fan(f- u)4 oe J2n(w— v)7 dydudw 


=-/- [ H(u,w) (S.C f ee fo) +S. — u =efopl Ol ener tee 
aan : (E.16) 


Carrying out the last integration with respect to w, results in 


Dina) =f Hluu tf) [Self ut fo) +Se(f-—4=fo)]du  (B.17) 


Evaluation of H(u,w) which is defined in Equation E.14, can be accomplished 
by performing a double Fourier transformation on n(r, 7). The single Fourier trans- 


form of n(r,+) is 


/ n(tachen 9" dn = / hy (y — a) / h,(r — a)e~??*"" dra, (a)da (£.18) 
Since 
/ hy(r — a)je~??""" dr = H, (u)e??*"* 


transforming again the left hand side of equation E.18, results in 


/ hr (y — aje??"" 7 dyda 


or 


HG) Gee [ Gg, (cr)en228 (4) do (E.19) 
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However, being a periodic deterministic signal, g,(a@) can be expressed in terms of 


its exponential Fourier series expasion as follows 


aa) = 3 Gyeren 


n=— co 


so that 


/ q.(a)e~2?7("-")*@da = c.f el2t(nRy—utw)a dy (E.20) 


— n=— oo 


Evaluation of the C’, coefficients for a particular case of interest is worked out in 


Appendix F. Combining Equations E.19 and E.20 yields 


H(u,w)= )> C,H, (u)H, (—w)6(nR, —ut+w) (Bro) 


n= -— 0 


From equation E-21, it is clear that 


H(u,u+v—f)= H,(u)H,(f —v—u) 5 Crd(nk, +u— f) (az) 


Equation E.17 can hence be further simplified to yield: 


Divx) =f . C,H, (u)H,(f —v—u)é(nRk, +v — f)- 


n=-— co 


(So(f + fou) +S.(f — fo —u)|du 


a Cn(nk, ty — 7) f ea oe 


n=- 60O 


(E.23) 


[So(f + fo —u) + S.(f — fo — u)|du 
Observing that in the integrand of equation E.23, u is the variable of integration, 


we can define 
Hf) & [ Ee eee) ne) S. (f= fy — u)ldu (B24) 
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Substituting now D(r,~) in equation E.1, with the use of equations E.23 and E.24 


yields: 


a | / H,(f)Hi(-v)P.(-f)P.(¥) S> C.6(nR, + u — f) HE (f,v)dudf 
- (£.25) 
Integrating Eq. E.25 over the variable v, we have: 


oO; =f Ss Cay) ae — f)P.(-f)P.(f —k,)H; (f, f —nk,)af 


n=— oo 


(E.26) 


From equation E.24, H; (f, f —nR,) can be seen to become: 
Hi(f,f—nR,) = / H, (u)H,z (nR, —u)(S.(f+fo—u)+S.(f—fo—u)|du (2.27) 


so that finally 


oo 


= > C, / ” Hy(f)Hr(nR, — f)P.(-f)P.(f — 2B): 


n=-— oo 


i H, (u)H,(nR, —u)[S.(f + fo —u) + S.(f — fo — u)|dudf 
i (E.28) 


12 


APPENDIX F 


THE EXPONENTIAL FOURIER SERIES 
OF q,(t) AND AVERAGE POWER OF THE PULSED JAMMER. 


Given q,(t) as defined in Chapter II, part C and re-illustrated as shown in 
Figure F.1 since q, (t) is a periodic signal, it is possible to expand q, (t) in terms of 
an exponential Fourier series with coefficients C,,, where: 


, 1 eT, . | T, 
oo =| | Aceuarn' (he a+ | c7 ef 27 t/Te dt (F.1) 
0 pT 


Let t; =t— — in the first integral and let t, =t— ae. in the second integral, 


so that 
pat (1-»9)T,/2 ue | 
an — = | A? ef24 (tt 0 S)/Te ge taf o2 pi2mn(ta+ yr, ‘dt, 


T 
a (1-)T,/2 


(F.2) 


Using properties of even and odd functions it’is possible to simplify Eq. F.2 to 














yield 
a (1-p) 
Cr = Ea 2a" el2nne/2 [ COs 2nnty dt, a 90? ef2mn(lt+o)/2 | Cos 2nTte aie 
T, 0 q | 0 Ty 
(1-p)Tq 
Be 25 ate ae sin2n7t, ae J” lend acre poate) sin 2nat, /T, 7 
T, 2g eel i 217) Dae 
= A2ypinn sinn7p +O(1L— peinn(ite) sinna(1 — p) 
NT p nm(1 — p) 
(F.3) 
From equation F.3, we observe that if n = 0, then 
Cy = A’ p+C7(1—p) (F.4) 
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bIl2 


(})§b uopfyouNY Ty einsry 


Br Lid bi(d-t)— 55- 


eV 


(7) *6 


74 


In the case where n 4 0, C,, can be reformulated as 


ayes nT p + (-1)"C2(1 — p) sinna(1 — p) 


C, =e" 
. nap nm(1 — p) 


yO (F.5) 


However, expanding sinna(1 — p) and taking into account that sinna = 0 and 
cosnm = (—1)", the second term of equation F.5 becomes —p? C? ee, so that 


we finally obtain 


sin n7p 


C, = e'"*"p?(A? — C?) ——— 0 6 
enn? p ( ae ei: (F.6) 
Summarizing the results above, we have 
A? p+C?(1—p ifn =0 
Cr a 2(A42 — A hs sin nm p if 0 (F.7) 
p*( Je — if 


The average power of the pulsed noise jammer can be obtained directly from 
the coefficients C,,. Alternatively, recall the expression for R,(t,t2) from equation 


C.1 and let t; =t, =t so that 


eee elt) — 2COs 2m fglt —t) [ h,(t — a)h, (t — a)q, (a)da 


or 


Ae ee . h* (t — a)q, (a)da (F.8) 


oO 
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Assuming specifically that h,(t) = Be~ ?*U(t), R, (t) can be obtained as follows, 


R, (t) =2/ [Be ?*-YU(t — a)]°q, (a)da 


CO oo 
ical | en yy Cel 27 lt-2)/Ts dy where rz=t—a 
0 n=-— co 


e(2B +j2nn/T,)z | 


= 2B? en jannt/t, Sail 
2 ; oe . 


: el2nnt/T, 
= 2B SSS 
Ly 2B + j2nx/T, 


2 el2nnt/T, 


= 25°) = 7 d,C-s 


el2mnt/t, 


2B —j2nz/T, + 20 Cn “2B + j2nz/T, 


BR? ‘ RP el2nnt/T, 
— — +2 . 
2 a > ater "2B + j2nnT, soe 
Since the C, coefficients have been shown to be of the form 
Ci fared. "2 fae va) (F’.10) 
Equation F.9 becomes 


Fea 2 


eIh te pjannt/T, 


Co — 
—-+2 R A, C, op} 
2B a cf Sal >) 3B + j2nr/T, 


ei(nxp+2nnt/T,) 
opt? Dh A,C,(\ RO, 
(2B)? = (2nx/T,)? e tan~!nnx/BT, 








= 25s 





— tan? a.) 
(2B)? (2h i 2 





C, cos ( 
= 2 ie op 2D fal (A, C, p) 


(F.10) 
It is obvious that R, (t) is periodic in t with period T,, since the average of a cosine 


signal is zero, we therefore have 





or 





P; = R,(t) = B[A’p + C*(1 — p)| (F'.11) 


APPENDIX G 


VARIANCE o% FOR FAST JAMMER 


In the case of fast jamming, the variance o% of the jammer as demonstrated 


in Eq. 3.18 simplifies to the following expression 
3 =, f lar RNP | iu wPS(f—wdudf (G2) 
where 


Dei Selo fo) e(2 — fo) 


In order to simplify the notation, denote the Fourier Transform and the inverse 
Fourier Transform operation by F{(-)} and F~!{(-)} respectively and let * denote 


the convolution operation. Define 


oS PV(N} (6.2) 
where 
v(t) = [|e (WP Sif —u)du (6.3) 


so that o4 can be reformulated as follows 


oF =O, [ ” |H, (£)PIP.(APV (A) af 


co 
oo 


of im (neie(ne | v(t\e%*!*dt of 
[one - 


=o fo v(t) fo late (AyPLP. (Alero af 


oo 


=o fo oF {IH (/PIPUNP} 2 


ia 


Recall from Chap. III, Eq. 3.17 that: 
P.(f) = G(f — fo) + G(f + fo) 
G(f)= sine (fie, 2 aee 
where for reasonable values of fy and 7, 
IP.(f)|? = P.(f)P.(-f) © |G(f -— fo) |? + |G(F + fo) (G.5) 
or equivalently 
2 Ty : . Q a : . 2 
ee eee @ sine“ (f — fo)T, + (4) since’ (f + fo)T, (G.6) 


The inverse Fourier Transform of this function can be easily obtained from tables 


(see Ref. 9), namely 


Fy ea ete ine te A @ (C2 eas) (Ga 
_ b 


jt], 
: 0 if ne 


The same considerations apply for |H;(f)|? since H,;(f) is the bandpass filter trans- 


where 


fer function specified in terms of shifted lowpass equivalents. Since H,(f) is signif- 


icant for f = +f over a restricted frequency band, 


|H:(f)|? + |Hie(f — fo)|? + |Hie(—f — fo)l? (G.8) 


where the lowpass equivalent H,,(f) is given by 


Btn 5 (G9) 
so that 
\Hio(f)|? = ; (a) (G.10) 
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Consequently, 
FSA ey = sent Peder (G.11) 


so that evaluation of F~* {|H;(f)|?|P.(f)|?} can be accomplished by using the 
convolution theorem [Ref. 9], that is 
P-* (1H (NPIP.UNP) = Fo CP} * PIP) 


i} A (= ea Els eit 20H e~ olt-7l ars Be poe) dr 
= b 


bee rT . 
= oa em e Olt-tla (=) cos 27 fo (t — t)dr 
bit Bree z 
= 7 te : 
8 ts A T, g(t, r)dr 
(G.12) 
where 
g(t, 7) = cos 27 fot cos 27 fo (t — 7) 
(G.13) 


1 a 
ES [cos 2m fot + cos 2mfotcos4afor +sin2afytsin 4m for] 


According to the value of |t — r|, Eq. G.12 can be further expressed as follows: 


. FAAP |P(AP}Hn+h (G.14) 
where: 
oT, f° T 
= af Pe elt) A € g(t, r)dr (G.15) 
I= os e TN (=| g(t, 7) dr (G.16) 
Se T, 


We can now proceed to the specification of v(¢) which is obtained from 
wart [ln wP sir - wan} 
=F" {|i (f)|?}- Fo" {5.(f)} (G.17) 


t 
= Be Fltla (= | cos 27 fot 


c 
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since 


F-'{\H,(f)|?} = Be7?"" cos2xfit (G.18) 


Substituting Eq. G.4 with Eq. G.14 and G.17 yields 


_ t 
On a CoB | e Bitla (=| cos 2m fot (I + I, ) dt 


— oo Cc 


The limit of integration can be reduced to the interval [-T,,T,] since the function 


Anfyt(I, + I,)dt The limit of integration can be reduced to the interval [—T,,T,| 


Laie: 
since the function A (=) is zero except for |t| < T,, hence 
c 


a ' 
on. = BES | eo Bltly € cos 27 fot(I, + I,)dt (G.19) 
—t. 


e 


According to the value of t, two cases can be considered as follows: 
Case (1) T; Seba QO, 


In this case, we have 


bT O t 
L+h= on | ere) ( + =| g(t,r)dr + | emits’? (: — =) g(t,7)dr 
, —-Ts b Oo cub 


eres ae 
+ | Ce (: ~ x) t,T ar 
Teed 


Since Eq. G.13 shows that g(t,r) is a function of cos4afyr and sin4mfg7, it can be 


(G.20) 


seen that the integration of functions of the form e”” cos4afp7 and e””’ sin 47fo7T, 
as well as forms te”™’ cos4mfor and te™’ sin4af ot will yield terms involving the 
factors (m? + 167? f2) and (m? + 167? f2)? in the denominator of the result while 
the numerators remain bounded. The numerators are proportional to m, which in 


this case equals 6 or B. Since in a practical system, 


fo>d and, —— =. 
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the contribution of the corresponding terms are therefore not significant, so that 


Eq G.20 can be simplified as follows: 


bT, 1 . T 
I PR ain 2 t — bt br 1 oe d 
IL+he 8 3 008 T fo f i3 ( +2) T+ 


t Tp 
cs | emer (: = =| dr +e” vs a = =| ar (G.21) 
0 T, t Ty, 


Eq. G.21 can be evaluated by using standard mathematics tables |Ref. 10] and 


therefore 


i eet e~ bT, 


bT;, 1 
+h = 60s 27 fot ; = br = BT, + on T. (e7 bt a e)| : (Ge22) 
b b 


Case (2) 0 >t > —-T,, 


In this case we have 


b t 0 | 
IL+l,= = if ents” € + =) g(t, 7) + / ge? | vat € + =| g(t, 7)dr 
-T;, b t b 


+ ‘ eeu € = =) a(tsr)ar 
(G.23) 


By considerations similar to those presented in Case (1), Eq. G.23 can be approx- 


imated by 


~ Oto = r as r eh 
ILr+he 3 5 cos 2 fot e js (+z) arte e€ ea T 
aad 


so that 


bT,, 1 t et ents 
i Iw — — + —— id a G.25 
1 + te 3 cos 27 fot ; a oT, PT i 20° T, (e aiane ) ( ) 
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We can now proceed to the evaluation of 04 expressed by Eq. G.19 according 


to the range of t. Since T, < TJ), 


BC, 67, ° t il t et € : 
2 0 b Bt | Z 9 | — bt ! bt 


, BC, oT, [ ~Bt t 2 1 t eu | Cia ae eet 
8 ; : z My ae 7 Jo 6 OT, 677, 2677, ¢ : 


(G.26) 














and using cos’ 27 fot = > (cos 4n fot+ 1) the reasoning that gave rise to the simpli- 
fications in evaluating Eq. G.20 can be used here to neglect the terms involving 
(m? + 16? f2) or (m? + 167? f2)? in the denominator. Thus Eq. G.26 can be 
further simplified as follows: 
BOOT. ff t\ fr tet enim 
2 = Bt 1 eon cal a ome =, is — bt =e bt dt 
Z 16 oe Tr) lb" bh eT, 267, ee 


ae t 1 t et e7 OTs 
—-Bt 2 pat. - Oe ee) a an — dt bt dt 
+ | : @ =| ; TT RR twp te | 
(G.27) 








and using standard mathematics tables (Ref. 10) results in 


BoCooT, {131 if i 2 
, o CoOL, -~BT , -BT 


1 1 — 6.7 2-047. ~BT.,5~oT,; 
-g(itee™) are t EE 








B? «BT, |2(B—b)?T, | 2(B+6)?T, 
(G.28) 
* By ea ee eS 2(B? + 67) 
B?—6? (B? —b2)?T, b?T, |B2—62 (B? —6?)?T, 
e~ (8 ->)T-. e (8 +6)T, 
+o BpPT STE 
which by regrouping terms, obtains the final form 
Cod 2 1 L+e7 27s 
2 O0-+b -BT 
SS ee ————— 
7a i (aE ar) (1 8") — BR 
B? B(B? +8?) 
2c? j) | eee G.29 
vue 1) EaG — 2) bT,T,(B? — 82)? oar 
Be7~2T< e7 T- Be-#8T- eof 
+ OOT,T. (B +6)? Y O67 T, (B — 6)? |} 
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APPENDIX H 
STATISTICS OF Y; AND Mi, 


We first recall the definitions of N,, and Y, from Chapter IV as follows 


Y,; = [ y(t)c(t) * h, (t) cos 2nf,tdt (7.1) 
Ty 


Nien = | Nea(t)c(t) * Ar (t) cos 2x f,tdt (H.2) 


0 


Observe that Y, is a zero mean random variable since the jammer and the spreading 
code c(t) are assumed uncorrelated and furthermore, c(t) has zero mean. The 


variance o% of Y; has been developed in Chap. 3 and can be expressed as (see Eq. 











3.19) 
Cal, a 2 at: oe 
= oh 1+ 7) some ~sx)- BT, 
. B? B(B? + 0?) 


—O9 ee ee Ne 
ma 4 ae —b?) 67, T,(B? — &)? ead 


Ben = tser he Be-28T- ebT- 
2b7,T.(B +6)? | 26T,T.(B — 6) | 
Replacing now B and 6 in Eq. H.3 by 27R, and 27 RP, respectively thereby setting 


-- 


the 3 dB cutoff of the filters appropriately, with 





Cn z 
defined as the processing gain, a final expression for of is obtained as follows 
gi “2 J(G.) (H.4) 
where 
l—e *")\(l-— az e727" = G?(2e7 7" —1 1 
ea Eg la 


: (e5)) 





—-+ 


Z G?4+1 20 7% eae eos 
21(G?2 — 1)? 4x |(G.+1)? (G, —1)? 
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The determination of the variance of, of the random variable N,, can be ac- 


complished in a manner very similar to that used in the determination of 04. Ob- 


serving however that replacing in the jammer model the deterministic priodic signal 


q(t) with a constant unit amplitude function, namely q(t) = 1, setting Ry (t,,t2) 


equal to “*6(t2,t,) (see Eqs. 2.16 and 2.17) results in a simple AWGN noise model 


having PSD level N, /2. Thus all evaluations involving the jammer can now be used 


with the stated modificiations in order to characterize the noise contributions to 


the performance of the system. Therefore, (see Eq. E.13) 


n(7,7) = [ h,(r — ajh,(y — a)da 


and with 
h,(t) = Be 7*U(t) 
then 
° M in(7r,7) 
n(T,7) = Breer 8 | e7 7 * da 
a B -Blr-1 
2 


Thus, the auto correlation function R, (t, rT) of the noise Ni, (t) is 
No po-Bir—7 
aca “Hee cos 2rf.(r — 4) 
From the methodology of Chap. 3, Eq. E.14 becomes 
= - B -Blr-—y|,-s2rru ,j2r7w 
Pi Noze Oe ae eo” Samay 


Evaluation of Eq. H.8 can be shown to yield 


B 1 1 
Luo) 
(u, w) Slee tro] (u — w) 
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(H.5) 


(1.6) 


(H.7) 


Therefore, the variance o7, of N,, (see Eq. E.25) 


cine | , SHAS )H,(—v)P.(—f) P.(—v) 


‘iy H(u,utv—f)[(S.(f —u+f.)+S.(f —ut+ f.)|dudfdv 


(H.10) 
where 
B 2B 
Huuty— N= Met ge bau (utv—f) H.11 
-y2( 2), (H.11) 
~ 22 | BP + (2ru2)? —" 


Letting f=v—f, or f=v+ f, as appropriate in the two integrals, obtain 


=f No lH; (f+fo)P.(f+f0)I? Z( pare f N,|Hi(f—f.)P.(f—f) 2 Za 


(H.12) 
where | 
~ B 2B 
4(f) =f- 2 B+ (anus ef — u)duaf . (H.13) 
Z(f) can be evaluated further as follows 
_ oo = Bp f.-sie 
Af)= 9 B? (27 f)? *5-(f) = aF te Re(t) 
= . | eB) R (t)e~F2"! dt (H.14) 


oe Pr. t 
= 5| | P* cos an feat — | 8 Ecos 2a feat 
1) Oo c 


The two integrals in Eq. H.14 can be evaluated separately and with BT, = 27, we 


have 
1 
A(f) = 5-4'(f) 
where 
A a (Ca cos 27 fT, — 1) (2) — 1) — 2FT, sin 2nT. )| 


(H.15) 
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From Figure H.1, it can be observed that if R, << R,, then Z'(f) = 1 for f within 


the interval [—R,, R,]. In fact, Z(f) is significant only for |f| << R,. Hence, Eq. 


H.12 can be simplified as follows: 


2 No 


ones ie Hi(F + LPF bPa + f= LIP L)Pad] 


20 


= Me fle PF Pa 


From Eq. G.6, we have 
Pi Soe) - 16S FS) ear 


where 


Ci (4) sinc’ fT,. 


(H.16) 


(H.17) 


The term |G(f — 2f,)| has significant components outside the integration region 


[—R,, R,]. Furthermore from Eq. Gis 


PE a) lee) ce ee laa |Hz,(—f)/ 


and again, the term |H,,(f —2/f.,)| can be dropped from further consideration due 


to reasons similar to those discussed above in the context of |G(f — 2f,)|. Thus 


Me Le 
Ge, zs xe f \H,,(f)|? (2) sinc’ fT, df 
Re 
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Z'f) 





z 
S 
| 


Figure H.1 Function Z’(f) 
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and since |H,,(f)|? has 3 dB cutoff at f = +R,, the integral itself can be approx- 





imated by 
aa sinc’ (fT,)df ~ k = 0.903 
Tie) oe i lati. 
Finally 
NT, 
Or > k (7.18) 
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